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I. INTRODUCTION 



O 



> 
oo 
in 



Recently, there has been renewed interest in the pairing problem in neutron matter and neutron-rich nuclei. The 
superfluid properties of neutron matter are of importance in the study of neutron stars [1] , while pairing in neutron- 
■ rich systems is of relevance for the study of heavy nuclei close to the drip line [2] and the light halo nuclei [3] . Much 
Q\ \ effort has gone into calculating the superfluid energy gap in dilute neutron matter [4-9]. Most of these studies have 
0^ ' been carried out using pairing matrix elements given by the bare nucleon-nucleon (NN) interaction. Even though it 
is a long time since Clark et al. [9] showed that density and spin-density fluctuations must be included in the pairing 
interaction, and there has been much progress in that direction recently [10,11], we will in this contribution focus 
on selected properties of the pairing problem in infinite neutron and nuclear matter employing only the bare NN 
interaction. 

\Q • In this lowest-order approximation to the problem it has been found that all modern NN potentials give nearly 
identical results for the 1 S'o energy gap in dilute neutron matter. One aim of this work is to explain how this can be 
understood directly from the measured properties of the free nucleon-nucleon (NN) interaction. This is discussed in 
Section 2. 

In Section 3 we discuss, still employing the bare NN interaction, various properties of the pairing wave function. 
The pairing gap is determined by the attractive part of the NN interaction. In the 1 S'o channel the potential is 
attractive for momenta k < 1.74 fm^ 1 (or for interparticle distances r > 0.57 fm). However, the nuclear situation is 
somewhat different from that of the classical BCS with attractive potentials in the solid state, see e.g., the discussion 
t-H , in Ref. [12]. In the so-called weak coupling regime, where the interaction is weak and attractive, a gas of fermions 
may undergo a superconduncting (or superfluid) instability at low temperatures and a gas of Cooper pairs is formed. 
This gas of Cooper pairs will be surrounded by unpaired fermions and the typical coherence length is large compared 
. with the interparticle spacing, and the bound pairs overlap. The latter behavior defines also what we will mean with 
weak-coupling in this work. With weak-coupling we will mean a regime where the coherence length is larger than the 
interparticle spacing. In the so-called strong coupling limit, the formed bound pairs have only a small overlap, the 
coherence length is small, and the bound pairs can be treated as a gas of point bosons. One expects then the system 
to undergo a Bose-Einstein condensation into a single quantum state with total momentum k — [12]. For the So 
'' • channel in nuclear physics we may actually expect to have two weak-coupling limits, namely when the potential is 
. £h \ weak and attractive for large interparticle spacings and when the potential becomes repulsive at r — 0.57 fm. In 
these regimes, the potential has values of typically some few MeV. One may also loosely speak of a strong-coupling 
limit where the NN potential is large and attractive. This takes place where the NN potential reaches its maximum, 
with an absolute value of typically ~ 100 MeV, at roughly ~ l f m . These properties of the NN potential in the 1 So 
channel and their connection with the wave function of the paired state are discussed in Section 3. In that section 
we will argue, from the properties of the wave function and the calculated coherence length, that fermion pairs in the 
1 Sq wave in neutron and nuclear matter, will not undergo the above-mentioned Bose-Einstcin condensation, since, 
even though the NN potential is large and attractive for certain Fermi momenta, the coherence length will always be 
larger than the interparticle spacing. 

Concluding remarks and further perspectives for pairing in nuclear systems are given in Section 4. 
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II. PHASE-SHIFTS AND PAIRING GAP IN INFINITE MATTER 



The energy gap in infinite matter is obtained by solving the BCS equation for the gap function A(fc). 

where V(k, k') is the bare momentum-space NN interaction in the 1 5'o channel, and E(k) is the quasiparticle energy 
given by E(k) = \J (e(k) — e{k F )) 2 + A(fc) 2 , where e(k) is the single-particle energy of a neutron with momentum 
k, and k F is the Fermi momentum. Medium effects should be included in e(k), but we will use free single-particle 
energies e(fc) = k 2 /2m, where m is the neutron rest mass, to avoid unnecessary complications. The omission of such 
medium effects is also in line with our omission of screening contributions. Anyway, at the densities considered here, 
Brueckner-type calculations [7] indicate that in-medium single-particle energies do not differ much from the free ones. 
The energy gap is defined as A F = A(k F ). Eq. (f) can be solved by various techniques, some of which are described 
in Refs. [7,8]. In Fig. f we show the results for Ap obtained with the CD-Bonn potential (full line) [13], the Nijmegen 
I and Nimegen 11 potentials (dotted line and dashed line, respectively) [14]. The results are virtually identical, with 
the maximum value of the gap varying from 2.98 MeV for the Nijmegen I potential to 3.05 MeV for the Nijmegen II 
potential. The same insensitivity of the results to the choice of NN interaction was found in Refs. [4,7]. We will now 
discuss how these results can be understood from the properties of the NN interaction in the So channel. 

A characteristic feature of 1 Sq NN scattering is the large, negative scattering length, indicating the presence of of a 
nearly bound state at zero scattering energy. Near a bound state, where the NN T-matrix has a pole, it can be written 
in separable form, and this implies that the NN interaction itself to a good approximation is rank-one separable near 
this pole [15]. Thus at low energies we can write 

V(k,k') = Xv(k)v(k'), (2) 

where A is a constant. Then it is easily seen from Eq. (1) that the gap function can be written as A F v(k), where A_f 
is the energy gap. Inserting this form of A(fc) into Eq. (1) one easily obtains 

Numerically the integral on the right-hand side of this equation depends very weakly on the momentum structure 
of A(k), so in our calculations we could take A(k) « A F in E(k). Then Eq. (3) shows that the energy gap A^ is 
determined by the diagonal elements Xv 2 (k) of the NN interaction. The crucial point is that in scattering theory 
it can be shown that the inverse scattering problem, that is, the determination of a two-particle potential from the 
knowledge of the phase shifts at all energies, is exactly, and uniquely, solvable for rank-one separable potentials [16]. 
Following the notation of Ref. [15] we have 

A^(fc) = - fc2 + ^ sin f^ e -"W, (4) 
Ar k 

for an attractive potential with a bound state at energy E = —n 2 B . In our case kb = 0. Here S(k) is the 1 So phase 
shift as a function of momentum k, while a(k) is given by a principle value integral: 

where the phase shifts are extended to negative momenta through 5(— k) = —6(h). 

From this discussion we see that Xv 2 (k), and therefore also the energy gap A^, is completely determined by the 1 So 
phase shifts. However, there are two obvious limitations on the practical validity of this statement. First of all, the 
separable approximation can only be expected to be good at low energies, near the pole in the T-matrix. Secondly, 
we see from Eq. (5) that knowledge of the phase shifts S(k) at all energies is required. This is, of course, impossible, 
and most phase shift analyses stop at a laboratory energy E\ a b = 350 MeV. The 1 So phase shift changes sign from 
positive to negative at E^ w 248.5 MeV However, at low values of k F , knowledge of v(k) up to this value of k may 
actually be enough to determine the value of A^ , as the integrand in Eq. (3) is strongly peaked around k F . 

The input in our calculation is the 1 5'o phase shifts taken from the recent Nijmegen phase shift analysis [17]. We 
then evaluated Xv 2 (k) from Eqs. (4) and (5), using methods described in Ref. [18] to evaluate the principle value 



2 



integral in Eq. (5). Finally, we evaluated the energy gap Ap for various values of kp by solving Eq. (3), which is an 
algebraic equation due to the approximation A(fc) « Ap in the energy denominator. 
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FIG. 1. 'So energy gap in neutron matter with the CD-Bonn, Nijmegen I and Nijmegen II potentials. In addition we show 
the results obtained from phase shifts only, Eqs. (3)-(5), and the effective range approximation of Eq. (6). 



The resulting energy gap is plotted in Fig. 1 (full line). As the reader can sec, the agreement between the direct 
calculation from the phase shifts and the CD-Bonn and Nijmegen calculation of A^ is, to say the least, satisfying, 
even at densities as high as kp — 1.4 fm -1 . The energy gap is to a remarkable extent determined by the available 
1 S phase shifts. In the same figure we also report the results (dot-dashed line) obtained using the effective range 
approximation to the phase shifts: 

k cot S(k) = — - + ^-r k 2 , (6) 
a 2 

where ao = —18.8 ± 0.3 fm and ro = 2.75 ± 0.11 fm are the singlet neutron-neutron scattering length and effective 
range, respectively. In this case an analytic expression can be obtained for Xv 2 (k), as shown in Ref. [16]: 

^ k 2 + r -£{k 2 + a 2 ) 2 V K Pi 

with a 2 = —2/aro, f3\ w —0.0498 fm -1 , and fa ~ 0.777 fm -1 . The phase shifts using this approximation are positive 
at all energies, and this is reflected in Eq. (7) where \v 2 (k) is attractive for all k. From Fig. 1 we see that below 
k F = 0.5 fm" 1 the energy gap can with reasonable accuracy be calculated with the interaction obtained directly 
from the effective range approximation. One can therefore say that at densities below kp = 0.5 fm -1 , and at the 
crudest level of sophistication in many-body theory, the superfluid properties of neutron matter are determined by 
just two parameters, namely the free-space scattering length and effective range. At such densities, more complicated 
many-body terms are also less important. Also interesting is the fact that the phase shifts predict the position of the 
first zero of A(fc) in momentum space, since we see from Eq. (4) that A(k) — Apv(k) = first for S(k) = 0, which 
occurs at E\ a b ~ 248.5 MeV (pp scattering) corresponding to k « 1.74 fm -1 . This is in good agreement with the 
results of Khodel ct al. [8]. In Ref. [8] it is also shown that this first zero of the gap function determines the Fermi 
momentum at which Ap = 0. Our results therefore indicate that this Fermi momentum is in fact given by the energy 
at which the 1 5 , phase shifts become negative. 

Thus, the quantitative features of 1 So pairing in neutron matter can be obtained directly from the Sq phase shifts. 
This happens because the NN interaction is very nearly rank-one separable in this channel due to the presence of a 
bound state at zero energy. This explains why all bare NN interactions give nearly identical results for the 1 S'o energy 
gap in lowest-order BCS calculations. However, it should be mentioned that this agreement is not likely to survive in 
a more refined calculation, for instance if one includes the density and spin-density fluctuations in the effective pairing 
interaction like in e.g., Rcfs. [10,11]. Other partial waves will then be involved, and the simple arguments employed 
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here will, of course, no longer apply. Our reasoning here applies also only to a partial wave where the T-matrix 
(almost) has a pole, and we have neglected the fact that the phase shifts become negative at higher energies. 
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2. 1 So energy gap in nuclear matter calculated with the CD-Bonn potential compared with the direct calculation from 
np and pp phase shifts. Also shown are the results for neutron matter with the CD-Bonn potential. 



The calculation of the 1 Sq gap in symmetric nuclear matter is closely related to the one for neutron matter. In 
fact, with charge-independent forces, like the older Bonn potentials, and free single-particle energies one would, of 
course, obtain exactly the same results. However, the new potentials on the market are charge-dependent, in order to 
achieve high quality fits to both np and pp scattering data, and therefore we should in principle solve three coupled 
gap equations for neutron-neutron (nn), proton-proton (pp), and neutron-proton (np) pairing: 



i r 
T Jo 



dk'k' 2 V t (k,k') 



E(k>) ■ 



(8) 



where i=nn, pp and np, and the quasiparticle energy is still given by E(k) = y / (e(k) — e(kp)) 2 + A(fc) 2 , but the 
energy gap is now given by 



A(fc) 2 = A nn (fc) 2 + A pp (fc) 2 + A np (£;) 2 . 



(9) 



Solving these equations, both with the CD-Bonn potential and with the phase shift approximations we get the results 
shown in Fig. 2. For comparison we have in the same figure plotted the results for pure neutron matter with the 
CD-Bonn potential (dashed line) . From the figure it is clear that the phase shift approximation works well also in this 
case. As could be expected, the results are very close to those obtained earlier with charge-independent interactions 
[7]- 



III. FEATURES OF PAIRING CORRELATIONS IN INFINITE MATTER 



An important length scale of e.g., a neutron superfluid is the coherence length. From a microscopic point of view 
the coherence length represents the squared mean distance of two paired particles (a Cooper pair of neutrons) on 
top of the Fermi surface. The magnitude of this quantity affects several of the physical properties of a neutron star 
crust. First of all, neutrons paired in a singlet state form quantized vortices induced by the rotational state of the 
star. These can pin to the nuclei present in the crust, possibly leading to the observed sudden release of angular 
momentum known as pulsar glitches. The magnitude of the pinning force depends on the size of the vortex cores, 
which is equal to the coherence length of the neutron superfluid. A second question is how properties of the neutron 
superfluid change due to the inhomogeneous environment of a neutron star crust, a problem related to the average 
thermodynamical properties of neutron matter [1]. The typical dimension of nuclei in the inner crust of a neutron 
star is Rn ~ 4 — 6 fm. This number is, in an appropriate range of densities, comparable to the coherence length £ as 
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estimated from existing BCS calculations. Clearly, the coherence length represents a critical parameter by which one 
can establish the behavior of an inhomogeneous superfluid. It sets the scale for the possible spatial variation of the 
pairing properties of the system, and thus plays a role if some inhomogeneities are present in the system at a length 
scale comparable to it. 

The coherence length can be easily evaluated from the wavefunction 0(r) of the relative motion of the two neutrons 
in a Cooper pair, r being the relative coordinate of the two particles. The coherence length £ is given by 

2 j^r|<Kr)|V _ J~dkk*\dx(k)/dk\* 

* Jd3 r | 0(r) |2 / °° d fcfc2| X (fc)| 2 ' ^ > 

with x(k) being the wavefunction of a Cooper pair in momentum space. This equation is particularly suited for 
numerical computation, since the BCS equations for a uniform system are solved in momentum space, as discussed in 
the previous section. The wavefunction of the Cooper pair in momentum space is given by (apart from an unimportant 
normalization constant) 

xW - §§, (ID 

where A(fc) is the fc-dependent pairing gap, while E{k) is the energy denominator in Eq. (1). 

As mentioned in the introduction, the NN potential in the Sq channel yields actually two so-called weak-coupling 
limits. This happens when the potential is weak and attractive for large interparticle spacings and when the potential 
becomes repulsive at r = 0.57 fm. This corresponds to the Fermi momentum kp — 1.74 fm -1 discussed in the 
previous section in connection with the phase shift analyses. The other weak-coupling limit, i.e., when r is large and 
the potential tends to zero, corresponds to small values of the Fermi momentum. The pairing gap decays exponentially 
to zero in the low-density limit [8]. What can be thought of as a strong coupling limit takes place where the NN 
potential reaches its maximum, at roughly ~ 1 fm, see e.g., Ref. [15] for a discussion of various features of the NN 
potential. The pairing gaps in Figs. 1 and 2 have their maxima at the density which corresponds roughly to the 
maximum of the NN interaction in the 1 So channel. 

In the weak-coupling limits, we have that the bound pairs overlap in r-space, or stated differently, that the wave 
function in Eq. (11) is strongly peaked in momentum space at the value of the corresponding Fermi surface. The 
coherence length £ is in this case much larger than the typical interparticle spacing, see e.g., the discussion in Ref. [12]. 
As stated in the introduction, we will with weak-coupling mean that the gas of Cooper pairs will be surrounded by 
unpaired fermions and the typical coherence length is large compared with the interparticle spacing, and the bound 
pairs overlap. 

If we have a strong coupling limit, the Cooper pairs at the Fermi surface have only small overlaps, a fact which 
means in turn that the pair wave function in Eq. (11) extends further out in fc-space, or that the Cooper pair is more 
localized in r-space. The coherence length should then be small, of the order or smaller than the interparticle spacing. 

The wave function of Eq. (11) for the various possible coupling regimes is shown in Fig. 3. There we plot x(fc) for 
five values of kp, 0.03, 04, 0.8, 1.2 and 1.4 fm -1 , employing the pairing gap from the previous section obtained with 
the CD-Bonn potential. The phase-shift approximation or the Nijmegen potentials yield essentially the same results. 

Clearly, at low values, kp — 0.03 fm -1 in Fig. 3, of the Fermi momentum, corresponding to one of the weak-coupling 
regimes, the wave function is strongly peaked in momentum space. Similarly, for densities where the NN interaction 
changes from being attractive to repulsive, we have the other weak-coupling regime. The qualitative form of the wave 
function at kp = 1.4 fm -1 resembles much that at low densities. For densities corresponding to r-values where the 
potential is close to its maximum, kp ~ 0.7 — 1.2 fm -1 , one sees that the wave function in fc-space is much more 
spread out, possibly implying that the coherence length is smaller and that the Cooper pairs have only small overlaps. 

The fact that the NN interaction in the 1 S'o channel is large and attractive at certain values of kp (up to five times 
larger than the Fermi energy) and that the wave functions in Fig. 3 extend over several values of k, may lead one to 
conclude that one could speak of bound fermion pairs which can be treated as a gas of point bosons. One expects 
then the system to undergo a Bosc-Einstcin condensation into a single quantum state with total momentum k — 0, as 
discussed in depth in Ref. [12]. However, such a conclusion for singlet pairing in neutron or nuclear matter is wrong. If 
one calculates the coherence length using Eq. (10) for the above Fermi momenta, one finds that for all Fermi momenta 
the coherence length is much larger than the typical interparticle spacing. Eq. (10) gives £ = 388.0, 4.8, 5.2, 13.2 and 
53.5 fm for kp = 0.03, 04, 0.8, 1.2 and 1.4 fm -1 , respectively. Even the smallest values are of the size of the radius 
of nuclei found in the crust of a neutron star. If one also observes that screening effects yield even larger coherence 
lengths, see e.g., [19], one can conclude that for singlet pairing in neutron or nuclear matter, the gas of Cooper pairs 
has a typical coherence length which is large compared with the interparticle spacing, the Cooper pairs overlap and 
fermion exchange may become dominant. 
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From Fig. 3 one also notices that for the chosen Fermi momenta, the pair wave function does not vanish at k = 0. 
This behavior is easy to understand if we again employ a rank-one separable interaction. Eq. (11) reads then 



X(k) 



v(k)A(k F ) 



which at k = simplifies to 



X(0) 



^/(e(k)-e(k F W+v(k)^A(k F )^' 
v(0)A(k F ) 



Ve(M 2 +«(0) 2 A(fc F ) 2 ' 



(12) 



(13) 



where A of Eq. (2) is set equal to one. The v(0) part of the potential can in turn be determined directly from 
the scattering matrix at k = 0. In that limit the scattering matrix equals — ao, where ao is the scattering length, 
a = —18.8 ± 0.3 fm for the neutron-neutron potential. For a rank-one separable potential, the on-shell scattering 
matrix at k — is given by [15] 



T(k = 0) = -a = 



v(0f 



1 + 



(14) 



If the 1 So channel really has a bound state at k = 0, the denominator should diverge, which in turn means that the 
scattering length should be a = — oo. The fact that the scattering length is finite implies that v(0) is finite. Eq. (13) 
can be rewritten as 



x(o)MM s 



v(0f. 



A(fc F ) 2 (l- X (0) 2 
When k F — > 0, the gap behaves asymptotically as [8] 

A(k F ) ~ 8e(fc F )e- 1/7 " 2 



(15) 



(16) 



where 7 = —2kfa /n. Inserting Eq. (16) into Eq. (15) then yields that the wave function x(0) — > when k F — > 0. In 
the other weak-coupling limit when k F = 1.74 fm -1 , i.e., where the potential changes sign, v(k F ) = 0, which inserted 
in Eq. (13) shows that the wave function goes to %(0) — > when k F — > 1.74 fm -1 . 
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FIG. 3. Wavefunctions in momentum space, x(k) for five different values of the Fermi momentum, kF = 0.03 fm -1 , kF — 0.4 
fm -1 , k F = 0.8 fm -1 , k F = 1.2 fm" 1 and fc F = 1.4 fm -1 . 

Finally, at k = kp, one sees from Eq. (12) that x(k — &f) = 1, as also seen in Fig. 3. 
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IV. CONCLUSIONS 



In summary, we have shown that in infinite neutron and nuclear matter, owing to the near rank-one separability of 
the NN interaction in the 1 So partial wave, we are able to compute the 1 5o pairing gap directly from the NN phase 
shifts. This explains why all NN potentials which fit the scattering data result in almost identical 1 5'o pairing gaps. 
Our findings conform with the conclusions of Khodel et al. [8] and Carlson et al. [20]: The virtual bound state in 
NN scattering determines the features of nucleon pairing in that partial wave. Even though this result is not 
likely to survive in a more refined calculation, for instance if one includes polarization effects in the effective pairing 
interaction like in e.g., Rcfs. [10,11], one can argue that our results demonstrate that upper limits for the value of the 
energy gap and for the density where a 1 S'o neutron/nucleon superfluid can exist, can be set directly from the 1 So 
phase shifts, since the polarization term serves to cut down the value of the gap, and leave the upper density for this 
superfluid more or less unchanged. These polarization terms will also enhance the already large coherence lengths 
for the singlet pairs, and one can conclude that for singlet pairing in neutron or nuclear matter, the gas of Cooper 
pairs will be surrounded by unpaired fermions, the typical coherence length is large compared with the interparticle 
spacing, and the bound pairs overlap. 

Finally, we note that fact that a bound state or a virtual bound state can be used to determine the properties of 
pairing in a physical system, may be of use in studies of superfluidity and superconductivity in atomic gases, such as 
a spin-polarized 6 Li gas, recently studied by Stoof et al. in [21]. The scattering length of lithium is large and negative, 
as is the case for the 1 Sq state discussed here. Since this is a very dilute system one can then even use an effective 
range approach to the inter-particle interaction and determine the gap uniquely for such dilute systems, by simply 
employing a separable interaction of the form shown in Eq. (7) and discussed in Fig. 1. 

We acknowledge several discussions with Marcello Baldo, Brett Carlson, John Clark, Manuel de Llano and Eivind 
Osnes. 
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